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ABSTRACT 


A study was made to improve the design of rocket engine acoustic liners. Helmholtz 
resonators with variable cavity depths were tested with an unstable rocket combustor. 
Combustor pressure oscillations inside the resonators and in the combustor and the phase 
angles between them were measured. From these measurements resonator aperture 
effective lengths were calculated. A new correlation for the effective lengths at high 
mean gas velocities is presented. A nonlinear resonator model is evaluated and used to 
calculate the aperture resistance. The calculated resistances were correlated by an 
empirical function of the mean and oscillatory flows. 
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SUMMARY 

A study was made to improve the design of acoustic liners for rocket engines. 
Helmholtz resonator aperture effective lengths were determined from rocket combustor 
measurements and a nonlinear resonator model was derived, evaluated, and used to cal- 
culate aperture resistance. To facilitate combustor measurements, five Helmholtz reso- 
nators with remotely variable cavity depths and known cavity gas properties were tested 
with a hydrogen-oxygen rocket engine. The engine was operated near the hydrogen tem- 
perature screech stability limit to produce measurable small pretransition oscillations. 

A measurement of the phase relation between the oscillatory pressure on both sides of 
the resonator aperture gave the timed cavity depths for the resonators. Aperture effec- 
tive lengths were calculated from the timed cavity depths and were compared with the 
mean gas velocity through the aperture. A new correlation for the effective length cor- 
rection at high mean gas velocities is presented. 

The measured ratio of cavity to chamber oscillatory pressure was compared with the 
results of the theoretical nonlinear resonator model. The results compared favorably 
with experimental data, and the theoretical model was then used to calculate the resona- 
tors' resistance as a function of mean and oscillatory aperture flows. These theoretically 
calculated resistances were correlated by an empirical function of the mean and oscilla- 
tory flows. 


INTRODUCTION 

The use of acoustic damping to suppress rocket engine combustion instability has 
been extensively investigated and reported in the literature (refs. 1 and 2). The present 



design of the acoustic dampers is based on maximizing the acoustic absorption of an array 
of Helmholtz resonators mounted on the walls of the combustion chamber. One method 
for calculating the damping has been reported in reference 1. This method depends on a 
knowledge of the correct values for resonator acoustic resistance 6 and resonator aper- 
ture effective length in a rocket combustor. Both the 6 and l are functions of 

the cavity gas properties, the acoustic wave amplitude, the resonator geometry, and any 
mean gas flow, whether past or through the aperture. Up to the present time, the calcu- 
lation of 6 and l has depended on estimating wave amplitudes, gas flows, and gas 
properties in the combustor and using these values in empirical correlations which are 
based on results with cold air flow (refs. 3 to 6). 

Errors in estimating these values can lead to large errors in calculating the acoustic 
absorption. In an effort to eliminate some of the design errors, the objectives of the 
present experiment were (1) to determine the aperture effective length Z eff from phase 
measurements made in a resonator mounted in a rocket combustor and (2) to develop a 
nonlinear resonator model, similar to that presented in reference 7, verify the model 
with resonator pressure measurements, and use the model to calculate 6, the resistance. 

To obtain the phase measurements, five Helmholtz resonators with variable cavity 
depths were mounted on a rocket combustor. The resonators were continuously purged 
with helium to assure known cavity gas properties. The combustor was operated near 
the hydrogen temperature screech stability limit to produce measurable combustor pres- 
sure oscillations. By measuring the phase relation between the oscillatory pressures on 
both sides of the resonator aperture as the cavity depths were varied, the cavity depths 
corresponding to tuning (oscillatory combustor pressure frequency = Helmholtz resonance 
frequency) were obtained. In the same manner as presented in reference 8, the Z eff 
was obtained from the known values of tuned cavity depth, oscillatory combustor pressure 
frequency, resonator geometry, and cavity gas sonic velocity. 

The nonlinear resonator model is that presented in reference 7, modified for the 
presence of the helium purge flow. The model was used to calculate values of cavity 
pressure oscillation amplitudes corresponding to a range of values of combustor pressure 
oscillation amplitudes. To verify the model, values of cavity and combustor pressure 
oscillation amplitudes were measured and compared with the calculated values. The 
model was then used to calculate values of the ratio of combustor pressure oscillation 
amplitude to the component of the aperture velocity which is in phase with the pressure. 
The normalized value of the ratio is 6, the aperture resistance. 

The investigation reported herein was conducted at the Rocket Engine Test Facility 

of Lewis Research Center using a liquid hydrogen-liquid oxygen rocket combustor. The 

combustor had an inside diameter of 10. 78 inches (27. 2 cm) and was operated at the 
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nominal values of 300-psia (2.07x10 -N/m abs) chamber pressure, 12 500-pound 
(55 600-N) sea-level thrust, and 95 percent C* combustion efficiency. 
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SYMBOLS 


2 2 

A aperture area, ft (m ) 

C adiabatic constant 

C* characteristic exhaust velocity, ft/sec (m/sec) 

c gas sonic velocity, ft/sec (m/sec) 

c heat capacity, ft/°R (m/K) 

d aperture diameter, ft (m) 

f frequency, Hz 

g c gravitational constant 

l aperture length, ft (m) 

l e ff aperture effective length, ft (m) 

M average mass per cross-sectional area, kg/ft (kg/m ) 

JK molecular weight, kg/kg-mole 

o 

P total pressure, psf (N/m ) 

p static pressure, psf (N/m ) 

R nonnormalized resistance, kg/(ft )(sec); kg/(m )(sec) 

r dimensionless weight flow 

T total temperature, °R (K) 

U aperture velocity, ft/sec (m/sec) 

AU mean-to-peak velocity in aperture, ft/sec (m/sec) 
mean velocity in aperture, ft/sec (m/sec) 
mean velocity past aperture, ft/sec (m/sec) 

V cavity volume, ft^ (m^) 

W defined by eq. (11) 

w purge weight flow, lb/sec (kg/sec) 

oi normal acoustic absorption coefficient 

5 2 (c n /c 

(3 dimensionless cavity geometry 

y ratio of heat capacities 
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6 aperture effective length correction, ft (m) 

6 acoustic resistance of resonator 

p gas density lb/ft ^ (kg/m'*) 

a array open area ratio 

r time, sec 

<f> acoustic resistance of resonator array 

co angular frequency, 1/sec 

co Q resonant frequency, 1/sec 

Subscripts: 

av average 

cavity resonator cavity 

max maximum 

0 reference condition 

1 chamber side 

II cavity side 

1 chamber side of aperature 

2 cavity side of aperture 

+ velocity into cavity 

velocity out of cavity 
Superscript: 

’ nondimensional 


APPARATUS 

Injector 

A 421- element flat-faced concentric-tube injector was used for the tests; the 
hydrogen- oxygen injection area ratio was 5. 18. The injector faceplate was copper with 
a hexagonal element pattern covering 85 percent of the face. A faceplate view of the in- 
jector is shown in figure 1. 
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Figure 1. - Injector face. 


Engine 

The heat- sink engine used had an inside diameter of 10. 78 inches (27. 2 cm) with a 

cylindrical section 10 inches (25.4 cm) long and a converging section 4 inches (10. 16 cm) 

long. The contraction ratio was 3.00. The engine was operated at a nominal mixture 

9 

ratio of 5 and a chamber pressure of 300 psia (2.07x10 N/m abs). 


Resonators 

Five resonators in a single axial row were tested simultaneously. The resonators 
were equipped with hydraulically actuated pistons for simultaneously varying the cavity 
depths during a run. The cavity depth was varied from ^ to 2^ inches (1. 27 to 6. 35 cm) 
in 2 seconds. The aperture thickness was 3/16 inch (0.476 cm); the aperture diameter 
was 1/4 inch (0.64 cm); and the piston diameter was 7/8 inch (2.22 cm). In addition to 
the instrumentation discussed in the next section, each resonator was equipped with a 
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Figure 2. - Injector resonator nozzle assembly. 


helium gas bleed line. A sketch of the resonators is shown in figures 2 and 3(a), and a 
photograph of the piston assembly is shown in figure 3(b). 


Facility 

The Rocket Engine Test Facility of the Lewis Research Center is a 50 000-pound- 
(224 000-N-) thrust sea-level test stand equipped with an exhaust gas scrubber and muf- 
fler. The facility utilizes a pressurized propellant system to deliver the propellants to 
the engine from the storage tanks. The facility is operated remotely from a control room 
located 2000 feet (610 m) away. A more complete description is given in reference 1. 


Instrumentation 

There were two distinct sets of instruments used for the tests, those for monitoring 
engine operation and those associated with the resonators. The instruments for engine 
operation are the same as those discussed in reference 1. 

The resonator instrumentation included a Chromel-Alumel open ball thermocouple 
within the cavity for measuring cavity gas temperature. In addition, a 1/16- inch 
(0. 159- cm) inside- diameter line was connected to a gas bottle to sample the cavity gas. 
Two high-frequency pressure transducers were used with each resonator cavity. One 
was mounted in the cavity piston, as shown in figure 3(b), and the other was mounted on 
the inside of the combustion chamber and 2-k inches (6. 35 cm) circumferentially from the 
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PreSSur ^s^ High-frequency pressure transducer 


•3/16 in. (0.476 cm) 


7/8 in. 

(2. 22 cm) 


(a) Cavity piston. 



C-67-2096 

(b) Hydraulic ram with cavity pistons. Each piston contains a high- 
frequency transducer, a thermocouple, and a gas bleed line. 


Figure 3. - Detail of cavity piston and hydraulic ram. 



aperture (the closest it could be placed). In addition, a high-frequency transducer was 
mounted on the inside of the chamber at an axial position corresponding to a resonator 
but displaced 135°. This was done to determine whether a resonant chamber mode was 
standing or spinning. The high-frequency transducers had nominal resonant frequencies 
of 10 000 hertz and a frequency response that was flat to 6000 hertz as mounted. The out- 
put from the high-frequency channels was recorded on analog tape and was later analyzed 
using a series of band pass filters and a phase meter to determine the amplitude fre- 
quency distribution of both the cavity and combustion pressure oscillations as well as the 
phase relation between them. 


THEORY 

The theory section is divided into two parts; in the first part, l or aperture ef- 
fective length is calculated, and in the second, the nonlinear resonator model is derived. 


Determination of 

The Helmholtz resonance frequency of a cavity- orifice combination is the following 
(ref. 1): 


w o = c cavity^^ (1 

As shown in reference 8, when the Helmholtz resonant frequency is equal to the im- 
pressed wave frequency, the phase angle between the impressed oscillatory pressure and 
the cavity oscillatory pressure is 90°. If the cavity volume V is altered so as to vary 
<x) Q until the phase angle between cavity and impressed wave pressure is 90°, then the 
impressed wave frequency u> equals cu q . Since cv is known from spectral analysis of 
the combustor pressure oscillations, c cayity is defined by gas purging, and A is fixed, 
then, at the value of V corresponding to co = u> 0 , is defined. 

Nonlinear Resonator Model 

As an aid in formulating the model, a sketch of a single resonator with relevant 
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dimensions is shown in figure 4. The assumptions inherent in the model, as described in 
reference 7, are as follows: 

(1) Isentropic expansions and contractions 

(2) Perfect gas 

(3) Constant velocity along l and across d 

(4) y same in cavity and chamber 
The momentum equation for the aperture is 


M dU 

iI37 = Pl ' P2 


(2) 


where M is defined by the following: 


M = P a v l eii 
P 1 +p 2 




For isentropic compression, the following holds: 


-o : 


l/y 


Substituting equations (3), (4), and (5) into equation (2) gives 


(3) 

( 4 ) 


(5) 
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dU 

dr 


eff 

2g, 



p l~ p 2 

1/y 


' p 2 


\l/r 


'II, 


Mass conservation for the cavity is satisfied by 


dr 


= w + 


p^V 


( 6 ) 


(7) 


where w is the helium purge weight flow and U is considered positive when directed 
into the cavity. The addition of the w then is the only modification to the model pre- 
sented in reference 7. 

Substituting equation (5) into equation (7) yields 


_n: 

nr 

dT 


\l/y 


= w + 


' P 2 


yl/r 


AU 


(8) 


When U is positive with flow into the cavity, the C of equation (8) is Cp whereas for 

U_, c = c n . 

The loss in the resonator is ascribed to loss of velocity head of the exiting jet. 
Therefore, for U + , the following hold: 


>i = p t(‘ - w ?) 


2 \y/(y-i) 


O) 


P 2 



where W + is defined by 


W 


+ 


U 


+ 


V 2c P g c T i 


( 10 ) 


( 11 ) 


For U_ , the following is assumed: 


p l = P I 


( 12 ) 
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( 13 ) 


/ o\r/(r-l) 

p 2 = p ii( 1 - w ) 


where W is defined by 


W = 


U 


/ 2 V c T n 


(14) 


Substituting equations (9) to (14) into equations (6) and (8) results in the following, 

/ 2 \r/(y-i) 

p i( l -_ w !) 


for U + : 


dt^ 

dr 


_eff 

2g r 


t 


, i( 1 - w ?) 


2\v/(r-l) 


1 /r 


(15) 


and, for U : 


= yw 4 /y + M I/V p 

dT V P^/vM V \ C I 


II 


w = 

+ 


u 


2y 


v l/2 


,V 


dU 


l) ( P l) (r " 1)/2rC I 1/2r 


„ ( 2\ y ^ y ~^ 
p i-p n (i-w 2 _j 


dT 

Z eff 

/P \ 1//y 
( l\ + 

/ 2V' // ^ -1 ^ 

Pul 1 - w :j 

l/y 

> 

2 *c 

\ c n / 

C II 

r 



- 



(16) 


(17) 


(18) 
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dp n_yw c n /y /. „,2V/(v-i) 


dr V (l/y)-l V 

P II 


’n (‘ - w ?)’ 


(19) 


W =■ 


u 


,1/2 


2y Y * l)/2y c l/2r 


( 20 ) 


To normalize these equations, the following parameters are chosen: 


re 


eff 


Al 


p i = — 

P o 


0 = 


P 2=~ 

P o 


eff 


a 2 = ' Cn 


,i/y 


u' = -5. 


W 

A 






O c 
p o o 


J 


Substituting the preceding parameters into equations (15) to (20) gives, 


for \J' + : 


(21) 


dU 


+ 2 


dr’ y 


/ 2 \y/(y-l) 

PU 1 -^ ) I p 2 

/ 2 \y/(y-l) 

"x (> - w !) 


+ ( P 2) Vr 


(22) 


dP 


dr 


2 = yjS^) 1_(1 / y W 2 + y^U; 


(23) 


U* 


w = 

+ 


,i/y 


(y - i 


( p i) 


t\(y-l)/2y 


(24) 
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for U' : 


du : _ 2a 2 P 1 - : 

/ 2 \r/(r-l) 

PH 1 - W -J 


dr’ y 

( p i) i/y - 

/ o\r/(r-i) 

P 2 v - W -/ 

i/r 


= y/3ra 2 (p ^"D/r + y /3U’ (i _ W 2 )^ ^ 


( 25 ) 


(26) 


W 



As discussed in reference 7, a reasonable approximation for 


a 


is the following: 


(27) 


a 2 


/%V /2 „ ^o T n 

\ c i / *^n T o 


(28) 


Equations (22) to (27) were solved by a fourth- order Runge-Kutta method on a digital com- 
puter. The input pressure Pj was sinusoidal only. The normalized part of impedance 
of the resonator was obtained by dividing the maximum input pressure by the in-phase 
maximum oscillation U T so that the following equation resulted: 


AP i 
y AU 


in phase 


(29) 


PROCEDURE 

The resonators were tested in a rocket engine that was operating stably but near the 
instability boundary. In reference 1, it was shown that a hydrogen- oxygen rocket engine 
can be driven into full-scale instability or screech by sufficiently lowering the hydrogen 
inlet temperature. As the unstable transition temperature was approached from the 
stable regime, a low amplitude chamber pressure oscillation was found to be present 
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which remained until the abrupt screech transition point (ref. 9). It was under this pre- 
transition condition that the resonators were tested. 

The first step in the testing sequence was to evaluate the screech transition tempera- 
ture of the injector. The procedure for testing the resonator was to set an inlet hydrogen 
temperature about 15° R (9 K) above the stability transition temperature. The damping 
due to the five resonators was not sufficient to affect the stability limits. The helium gas 
bleed system was in operation at all times to fix the cavity gas composition. During 
engine startup, the helium purge was maintained at a high value to prevent the ignition 
spike from forcing hot combustion gases into the cavity. After steady- state combustion 
had been achieved, at some predetermined time in the run, the gas purge flow was re- 
duced to the desired level and the five hydraulically operated pistons were moved simul- 
taneously to increase the cavity volumes. During the piston motion, a valve in the gas 
sample system was opened and closed for a predetermined time and the gas sample stor- 
age bottles were filled to 300-psia (2. 07x10 -N/m abs) chamber pressure. The total 
run time was 2^ seconds. 


RESULTS AND DISCUSSION 

To determine the values of aperture effective length it was necessary to mea- 

sure the cavity and chamber oscillatory pressures, the phase relation between them, the 
cavity gas temperature, the cavity purge flow, and the cavity depth as a function of time 
for each run. 

Typical results for a run are shown in figure 5. The cavity depth was maintained 
constant as the mean chamber pressure rose to a steady- state value of 300 psi (2. 07x10 
N/m ). At r = 0.6 second, the hydraulically controlled pistons were withdrawn and the 
cavity depth was linearly increased. The corresponding cavity gas temperature also 
changed during this period as shown in figure 5(b). In the initial portion of the run, a high 
purge flow was maintained to prevent the ignition spike from forcing hot gases into the 
cavity. From r = 0. 2 to 0.6 second, the high purge pressure was gradually reduced to a 
value which was kept constant for the remainder of the run. The surge in temperature 
corresponded to the decrease in purge flow. At the screech transition point, the temper- 
ature of the cavity increased markedly. An example of the output of two high-frequency 
transducers (one in the resonator cavity, the other on the chamber wall), is shown in fig- 
ure 5(c). For this particular run, it was not possible to maintain the hydrogen injection 
temperature above the transition temperature and high amplitude instability resulted at 
t - 1.8 seconds. 

The phase relation between the two transducers is shown in figure 5(d). As the cav- 
ity depth was increased, the ratio of oscillatory chamber pressure frequency 
(5300 Hz) to cavity or Helmholtz resonant frequency was also increased and the phase re- 
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(c) High-frequency transducer signal. 



(d) Phase angle between cavity and chamber. 
Figure 5. - Typical results for a run. 


lation between the pressure oscillations changed as shown. Comparison of the outputs 
from the two transducers located 135° apart at the same axial position indicated that the 
chamber oscillation pattern or mode was standing rather than spinning. Thus, the phase 
angle between transducers located 2-^ inches (5. 35 cm) from the aperture and inside the 
cavity corresponded to the phase relation across the aperture with no error due to the 
2-^-inch (5. 35-cm) displacement introduced. At a phase lag of 90°, cavity reactance be- 
came zero and the wave frequency was equal to the cavity resonant frequency. This is 
the tuned condition. 

An important step in analyzing the data was to obtain spectral analyses of the output 
of the high-frequency transducers as a function of cavity depth since this gives the value 
of ye or u>. A typical spectral analysis of the output of a high-frequency transducer 
mounted on the combustion chamber is presented in figure 6. The pretransition period is 
characterized by a frequency of 5300 hertz. After the screech transition, the amplitude 
increased markedly and the frequency shifted to 5700 hertz, corresponding to the second 
tangential mode. 

15 
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2468 lOxlO 3 
Frequency, Hz 

Figure 6. - Spectral analysis of signal from chamber high- 
frequency transducer for pretransition operation. 

Calculation of Aperture Effective Length 

The aperture effective lengths of the resonators were calculated in the same manner 
as described in the section Determination of Gas properties and the mean aper- 

ture velocities were defined by the helium purge, and the aperture area was fixed so that 
the cavity volume and aperture effective length controlled the resonant frequency. As the 
cavity volume or depth was increased, the phase angle between the pressure oscillations 
increased, as shown in figure 5. At a phase angle of 90°, corresponding to a 1/4-cycle 
lag, f wave was equal to the % elmholtz (ref. 8). Since f waye was known (5300 Hz) and 
cavity volume V was known as a function of time, the aperture effective length could be 
calculated. 

Before presenting the results of the calculation, a brief discussion of the ef- 

. fects of mean flow through and past the aperture is in order. The effect of mean flow 
past an aperture on effective length was reported in reference 5 and is presented in fig- 
ure 7. A mean flow of 394 feet per second (120 m/sec) reduces the aperture mass or ef- 



20 40 60 80 _ 100 120 140 

Mean flow past aperture, U, m/sec 

Figure 7. - Effect of mean flow past on aperture mass or t eff . Aper- 
ture length l = 0.3 centimeter; aperture diameter d = 0.8 centi- 
meter. (See fig. 3 in ref. 5.) 
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Mean flowthrough aperture, U, cm/sec 

Figure 8. - Effect of aperture flow through on aperture 
reactance. Aperture length l = 0. 0508 centimeter; 
aperture diameter d = 0.357 centimeter. (See fig. 13 
in ref. 6.) 




0 20 40 60 80_ 100 120 140 


Mean flow past aperture, U, m/sec 

I I I I I J I I 

0 1 2 3 4 = 5 6 7 

Mean flow through aperture, U, m/sec 

Figure 9. - Comparing effects of mean flow past and mean flow through on aperture 
mass. 
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fective length by a factor of 0.625 or reduces it to a minimum or saturated value of 0. 375 
times its nonflow value. A separate study, reported in reference 6, indicated that mean 
flow through the aperture would also have a strong effect on Z e ^. The flow-through ef- 
fect is presented in figure 8. Changes in acoustic reactance for this experiment are 
equivalent to changes in effective length. The results indicate a similar reduction with 
the flow through to that with flow past. In order to compare the two effects, both results 
are plotted as the ratio of with flow to without flow and are shown in figure 9. 

The agreement between the two effects indicate a similarity in l ^ reduction. Compar- 
ison of the flow velocities required for the same effect indicates the flow through is ap- 
proximately 20 times as important as flow past. In the rocket engine environment dis- 
cussed in this report, the mean flow past was approximately 1000 feet per second (305 
m/sec) and the mean purge flow through the aperture was approximately 250 feet per sec- 
ond (76 m/sec). Either effect, or a combination of the two, was sufficient to reduce the 
aperture effective lengths to their minimum or saturated values. Therefore, the results 
presented in figure 10 are Z ^ values which are saturated and thus independent of flow. 



Figure 10. - Aperture effective length at large flow-past and 
flow-through velocities. Velocity past, 1000 feet per second 
(305 m/sec); velocity through, 250 feet per second (76 m/sec). 


The data are shown as points, and the dashed lines indicate values of l ^ predicted by 
empirical correlations shown in equations (30) to (34). 

The expression for l ^ without mean flow past or through the aperture is 


eff 


= l + 5 


(30) 


where 


6 = 0. 85 d 1 - 




cavity area 


(ref. 1) 


(31) 
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The topmost dashed line in figure 10 is based on equation (30). 

The expression generally used to calculate the saturated or minimum l is taken 
from figure 9; it is as follows: 

l eff' (saturated) = 0. 375 (l + 5) (32) 

The bottom dashed line in figure 10 represents this expression. 

Based on a limited quantity of data, the results of reference 8 indicate that the fol- 
lowing expression may be valid: 

Z e ff(s a turated) = l (33) 

The line representing equation (33) is indicated as such in figure 10. 

Of these three correlations, equation (33) seems to represent the data best; how- 
ever, it is still not satisfactory. 

Equation (32) was developed from experiments using small aperture thicknesses 
(Z = 0.010 to 0. 118 in. or 0. 05 to 0. 3 cm) so that the aperature effective length correc- 
tion factor 5 was significantly larger than the thickness Z. Therefore, the data cannot 
distinguish whether flow reduces Z + 5 or only 5. Reference 10 presents a theoretical 
argument showing that, for thin orifices (Z ~ 0), the aperture effective length is reduced 
by turbulence to 0. 375 times its original value. Consequently, a proposed expression for 
the saturated Z e ^ for thick- wall resonators is 

Z eff (saturated) = l + 0. 375 6 (34) 

The dashed line, representing equation (34) and indicated as such in figure 10, agrees 
best with the data, indicating that equation (34) should be valid for the high mean flows 
associated with rocket engine liners. 

Discussion of Nonlinear Resonator Model 

As indicated in the INTRODUCTION, an important step in the calculation of 6, the 
acoustic resistance, is the formulation and validation of a nonlinear resonator model. 

The model is formulated and presented in the section Nonlinear Resonator Model. A 
typical result of using the model is shown in figure 11. The output presented in the 
calculated cavity pressure as a function of normalized time. The input parameters are 
Pj max = r = 0, y = 1.2, j8 = 0. 10, and ~cft = 1. 00. These values were chosen to 

correspond to those presented in reference 7. The results reach a steady state after the 
first cycle at a peak-to-peak value of 0. 365 chamber pressure with a direct- current shift 


19 



Figure 11. - Output of nonlinear model cavity pressure response. Ratio of frequency to 
Helmholtz frequency f/fnelmholtz * 0- 855; ratio of heat capacities y - 1.2; ar • 1.0; 
dimensionless cavity geometry p * 0. 1; dimensionless weight flow r - 0; maximum 
chamber pressure P^ max = 0. 10. 


of 0. 012 chamber pressure. This shift is due to the assumptions inherent in the model 
and the nonlinearity of the model. 

Both the velocity and cavity pressure outputs as a function of normalized frequency 
are shown in figure 12. The value taken as representative of the output is one -half the 
peak-to-peak value. These representative values are indicated by AU' and APg. The 
normalized frequency is the ratio of wave frequency to Helmholtz resonator frequency 
u>/u> . For purposes of comparison, the results of the same configuration as presented 
in reference 7 are superimposed. The good agreement indicates the comparability of the 


_c\j 

Q_ 

<1 



Figure 12. - Normalizedoscillatory velocity 
and normalized oscillatory cavity pressure 
as function of ratio of wave frequency to 
Helmholtz frequency. cir= 1.0; dimension- 
less cavity geometry p = 0. 10; dimension- 
less weight flow r = 0; maximum chamber 
pressure PJ max = 0. 10. 
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Dimensionless 
weight flow, 



Ratio of wave frequency to Helmholtz 
frequency, colu) 0 


Figure 13. - Effect of mean flow through on 
resonator response as function of ratio of 
wave frequency to Helmholtz frequency. Ra- 
tio of heat capacities y ■ 1. 6; - 1. 0; 

dimensionless cavity geometry (3 = 0. 02; 
maximum chamber pressure max = 0.02. 
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Figure 14. - Comparison of theoretical and ex- 
perimental oscillatory cavity pressures as 
function of oscillatory chamber pressures 
for 90 c phase lag. Frequency f = 5300 hertz. 


integration procedures. The velocity maximizes in the region of co/cJq, whereas the 
cavity pressure maximizes at frequency ratios of less than 1. The results are similar to 
the response of a forced-motion mass-spring-damper system with the shift in the peak 
cavity pressure response a function of the resonator damping. 

The effect of purge flow on the oscillatory cavity pressure is presented in figure 13. 
The results are presented as a function of the dimensionless frequency cu/ cu for various 

o o 

values of dimensionless weight flow r defined as (w/A)/p o c q . The values of y, /3, a , 
and PI are representative of experimentally measured values. As r is increased, 

J. ^ XXI 

the cavity pressure response curve becomes lower and flatter in the same manner as 
would be expected for large increases in the aperture damping. 

In order to validate the resonator model, the cavity gas composition, temperature, 
cavity depth, and oscillatory combustion pressures were converted to input parameters 
and the oscillatory cavity pressure calculated. Presented in figure 14 is a comparison of 
the corresponding theoretical and experimental oscillatory cavity pressures. Because of 
the difficulty in measuring the oscillatory pressures at a specific time in the run, a phase 
angle of 90° between combustor and cavity oscillatory pressures was used to define the 
only time at which theory and experiment were compared. The correspondince between 
theoretical and experimental values is indicated by the dashed lines in the figure. 

Although the agreement was fairly good and tended to suggest the validity of the 
model, the lack of better agreement was likely because of the following factors: 

(1) There is no provision in the model for including the effects of mean flow past. 
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(2) Inherent in the model is the assumption that the y within the cavity is the same 
as in the combustor. This assumption is in error by 30 percent, which can introduce 
changes in predicted values by up to 30 percent, but no change in trends. 

(3) The validation of the model depends on correlating the temperatures, cavity 
depths, and oscillatory pressure amplitudes all at precisely the same time. This is 
particularly difficult with pressure oscillations that demonstrated variations over 30 per- 
cent during a short time period. 

(4) The inability to measure the oscillatory combustor pressure immediately adjacent 
to the resonator aperture. The relation between the oscillatory pressure at the resonator 
aperture and the closest measurement point (2-g in. or 6. 35 cm away) could not be pre- 
cisely calculated and, therefore, the two were assumed equal. 

Considering all of the possible sources of error, the agreement shown can be considered 
a limited verification of the model. 


Calculation of Resistance 0 

If it is assumed that the comparison between theory and experiment shown in fig- 
ure 15 represents a limited verification of the resonator model, it is then reasonable to 
use the model to compute values of the resistance for various conditions with and without 
mean through flow. 

The results of a number of variations in oscillatory chamber pressure and mean flow 
are presented in figure 15. The variations in oscillatory chamber pressure cause varia- 
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Figure 15. - Correlating theoretically calculated acoustic 
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tions in oscillatory flow velocity. The acoustic resistance, defined in the section Theory 
as the ratio of normalized oscillatory chamber pressure to normalized in-phase aperture 
velocity, is plotted as a function of an empirical correlation of mean and oscillatory 
flows. The empirical correlation is arrived at as follows: the results of the study re- 
ported in reference 11 indicated that the nonnormalized acoustic resistance R can be ex- 
pressed as 


R = 1 . 5 p Q U (35) 

when U is significantly larger than AU, the oscillatory velocity. When a pressure loss 

2 

term of the form p Q U is assumed, in the absence of mean flow, R can be expressed as 

R = P Q AU (36) 

An empirical combination of equations (35) and (36) is 

— = 1. 5 U + AU (37) 

^o 

Normalizing equation (37) gives the following: 

6 = 1,5 U' + AU’ ' (38a) 

Equation (38a) is presented as the upper dashed line in figure 15. The computed values 
lie somewhat lower as indicated by the symbols. The half-filled symbols, representing 
oscillatory flow only, seem to agree with the empirical dashed line better than the points 
representing both mean and oscillatory flows. As the flow velocities increase, the re- 
sults deviate increasingly from the empirical correlation because the pressure loss term 
used in the model is not p Q U but the more complicated expression of equation (9). The 
difference between these two pressure loss terms increases with increasing velocity. A 
more suitable expression for the resistance which better fits the results is 

0=0. 8(1. 5 U* + AU') (38b) 

This expression is shown as the lower dahsed line in figure 15. The use of the empirical 
model, as represented by equation (38b), would greatly simplify previous liner design cal- 
culations (e. g. , ref. 1). 
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Use of the Empirical Resistance Equation 


The use of the empirical resistance equation (eq. (38b)) can be demonstrated using a, 
the acoustic absorption coefficient. Under certain conditions, a can characterize the 
liner damping (ref. 12). In order to use a, the resistance of a single resonator 6 must 
be converted into the resistance of an array of resonators <f>. This is accomplished by 
dividing 6 by the liner open area ratio a so that the following equation results: 

4> = — (39) 

a 

The equation for a for a liner that is tuned is 

a = 4 — (40) 

(* + l) 2 


Dimensionless mean- 



Figure 16. - Effect of mean flow on acoustic absorption coefficient a. 
Liner is tuned; open area ratio a = 0. 10; acoustic resistance of 
resonator 0 = 0.8 (1.5 U' + AU 1 ). 


Combining equations (38a), (39), and (40), the variation of a as a function of U’ for two 
values of AU' is calculated and presented in figure 16. The mean flow can decrease the 
absorption coefficient by a factor of two, although its effect is lessened for higher values 
of oscillatory flow AU'. The results shown in figure 16 are important since they demon- 
strate that when a large value of gas purge flow is required to fix liner cavity proper- 
ties, it can seriously reduce liner absorption. 
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SUMMARY OF RESULTS 

Comparison of the calculated values of aperture effective length Z gff with empirical 
correlations indicate that the following expression should be used in a rocket engine liner 
with high mean flow past or through the aperture: 


Z e ff = l + 0. 375 6 



where l is aperture length, 6 is aperture effective length correction, d is aperture 
diameter, and A is aperture area. 

There is reasonable agreement between the theoretically calculated values of cavity 
oscillatory pressure and those which were experimentally measured. Using the theo- 
retical model, values of the acoustic resistance, 0, were calculated. Comparison of the 
calculated values of 0 with an empirical correlation indicated that the following expres- 
sion best represented the data: 

e= 0. 8(1. 5 U’ + AU’) 

where U' is nondimensional mean velocity in aperture and AU' is nondimensional mean- 
to-peak velocity in aperture. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 12, 1968, 

128-31-06-05-22. 
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